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THE STAGING DYNAMICS OF A PROPOSED 

SPACE SHUTTLE CONFIGURATION 

By David A. Hamilton 
Manned Spacecraft Center 


SUMMARY 


A mathematical model of the staging dynamics of a proposed space shuttle con- 
figuration was developed to evaluate a ram-activated linkage separation mechanism. 
Included in this model was the kinematics and dynamics of the staging mechanism, 
thrust forces, thrust vector control, rigid-body dynamics, and structural dynamics of 
both the booster and orbiter stages of the configuration. The model was incorporated 
into a computer program so that the staging maneuver could be simulated. The results 
of simulations of the staging maneuver demonstrated that the link- ram separation de- 
vice provides satisfactory separation at nominal staging conditions. 


INTRODUCTION 


In an attempt to economize the U. S. space program, a completely reusable earth 
orbital launch vehicle has been under consideration by the National Aeronautics and 
Space Administration (NASA). One concept of this vehicle is a two-stage configuration 
in which the upper stage, the orbiter, is mounted ’’piggyback" on the lower stage, the 
booster. Both the booster and orbiter are winged vehicles. The configuration makes 
a vertical, rocket-type ascent, and each vehicle makes a horizontal, airplane-type 
descent. The booster delivers the orbiter to a desired altitude, then separates, flies 
back to earth, and lands. The orbiter continues thrusting until earth orbit is achieved. 

The separation of the orbiter from the booster is a critical maneuver. The ma- 
neuver must occur safely, that is, the linear and angular motion of each vehicle must be 
such that no recontact between the vehicles occurs, and excessive heating because of 
engine plume impingement should not occur. A separation mechanism must be designed 
that will ensure safe separation. The loads generated on each vehicle by the mechanism 
must be within design limits, and the operation of the mechanism must be such that 
orbiter velocity losses during the maneuver are as small as possible. 

Many candidate separation mechanisms have been proposed by NASA and support 
contractors. Because the testing of all proposed mechanisms is impractical, they are 
evaluated primarily by analytical techniques. These techniques involve developing 
mathematical models of the staging system. The models include the equations of mo- 
tion of the booster and orbiter and the analytic description of the various forcing 


/ 



functions that affect the motion of each vehicle. The mathematical model is usually 
incorporated into a computer simulation so that the staging maneuver can be simulated. 

This thesis is concerned with the development of a mathematical model for pre- 
dicting the interaction dynamics of the booster and orbiter during the separation se- 
quence for a ram-activated linkage separation mechanism. 

The link- ram system was not selected as the prime system by any of the NASA 
support contractors; thus, only limited analyses were conducted to evaluate it. The 
development in this thesis differs from previous analyses in that it includes modeling 
of the full six-degree- of-freedom motion of each vehicle, the structural representa- 
tion of each link member, and the structural response of the booster and orbiter. 

The mathematical techniques utilized in modeling the separation dynamics of the 
"piggyback" shuttle configuration are also applicable to the parallel burn shuttle con- 
figurations being considered by NASA. 


SYMBOLS 


A 


BEV 

[BMG], [OMG] 


BRMV 1 


brmv 2 


CX, CY, CZ 
DLX. 

l 


DLY. 


dly ri 

DTBl , DTB2, DTB3 
DX, DY, DZ 
EBX, EBY, EBZ 
EBXI, EBYI, EBZI 
E 


link cross-sectional area 

vector location of booster engine 

matrix of natural frequencies of each vehicle 

force vector on booster produced by ram mechanism 1 

force vector on booster produced by ram mechanism 2 

effective velocity gains 

axial deflection of ith link member 

lateral deflection of ith link member 

lateral deflection of ram mechanism 1 

booster Euler angle errors 
gimbal command gains 
booster error signals 
booster error signal integrals 
modulus of elasticity of link members 
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ETBX, ETBY, ETBZ 

FI, F2, F7, F8 

FBE 

FLB. 

1 

FLKB 

FLKO 

FLO. 

l 

FRAM 

FRBj 

frb 2 

FRM 

FROj 

fro 2 

FRY 1 


booster coordinate system components of Euler angle errors 

forces acting at end points of link element 

external forces acting on booster 

force vector at booster ith link attach point from link 
deflections 

total force on booster from link forces 

total force on orbiter from link forces 

force vector at orbiter ith link attach point from link 
deflections 

magnitude of ram force 

lateral force vector on booster from ram mechanism 1 

lateral force vector on booster from ram mechanism 2 

maximum magnitude of ram force 

lateral force vector on orbiter from ram mechanism 1 

lateral force vector on orbiter from ram mechanism 2 

lateral force on orbiter from ram mechanism 1 


, ^ 

FSMB total force vector on booster 



FSMO total force vector on orbiter 

G shear modulus of elasticity of link element 


GBCX., GBCY., GBCZ. 

J J J 


gimbal angle commands of jth booster engine 


[GBE(I, J)]j 


transformation matrix from booster to jth booster engine 
coordinate system 


[GBEA] 


booster Euler rate transformation matrix 


[GBL(I, J)]- transformation matrix from booster to ith booster link coor- 

1 dinate system 

[GBOL(I, J)] . transformation matrix from ith booster link to ith orbiter 

1 link coordinate system 
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GBXD., GBYD., GBZD. 
J J J 

[GIB (I, J)] 

[GIO(I, J)] 

[GOE(I,J)j. 

[GOEA(I, J)] 
[GOL(I,J)]. 

I 

[INRB(I, J)] 

[lNRO(I, J)] 

IXXB, IYYB, IZZB 
IXYB, IXZB, IYZB 
IXXO, IYYO, IZZO 
IXYO, IXZO, IYZO 
ij,k 

J 

KRL 

KRR 

i 

[MBG], [MOG] 

MBT 

MOT 

M 1 

M 2 

OEV 

4 


gimbal rate commands of jth booster engine 

transformation matrix from inertial to booster coordinate 
system 

transformation matrix from inertial to orbiter coordinate 
system 

transformation matrix from orbiter to jth orbiter engine coor- 
dinate system 

orbiter Euler rate transformation matrix 

transformation matrix from orbiter to ith orbiter link coor- 
dinate system 

cross-sectional moment of inertia of link element 

inertia matrix of booster 

inertia matrix of orbiter 

booster moments of inertia 

booster products of inertia 

orbiter moments of inertia 

orbiter products of inertia 

unit vectors along X-, Y-, Z-axes 

polar moment of inertia of link element 

lateral stiffness of ram mechanism 1 

engine roll control constant 

length of link element 

generalized mass matrix of each vehicle 

total booster mass 

total orbiter mass 

number of booster engines 

number of orbiter engines 

vector location of orbiter engine 



ORMV 1 

ormv 2 

PBSR 1 


pbsr 2 

POSB. 

l 

POSO. 

POSR^ 

posr 2 

RX, RY, RZ 

SB 

SBD. 

SBD R 1 

SBD R2 


SBO 

SEPV. 

SO 

SOD. 

1 


SODrI 

SODR2 

SPXI. 

1 


SRMj 

SYI 1 


force on orbiter produced by ram mechanism 1 

force on orbiter produced by ram mechanism 2 

vector from booster center of gravity (c. g. ) to ram 1 attach 
point on booster 

vector from booster c. g. to ram 2 attach point on booster 

vector from booster c. g. to ith link attach point on booster 

vector from orbiter c. g. to ith attach point on orbiter 

vector from orbiter c.g. to ram 1 attach point on orbiter 

vector from orbiter c.g. to ram 2 attach point on orbiter 

angular rate gyro gains 

inertial position vector of booster c. g. 

booster structural deflection vector at ith link attach point 

booster structural deflection at ram 1 attach point 

booster structural deflection at ram 2 attach point 

vector from booster c.g. to orbiter c.g. 

vector along ith link from booster to orbiter attach point 

inertial position vector of orbiter c. g. 

orbiter structural deflection vector at ith link attach point 

orbiter structural deflection at ram 1 attach point 

orbiter structural deflection at ram 2 attach point 

initial axial length of ith link 

vector along ram 1 from booster to orbiter attach point 


initial lateral deflection of ram mechanism 1 



T4, T6, T10, T12 
TBRMj 

tbrm 2 

THMB 

THMO 


THRVB 


THRVO 

TLB. 

1 

TLKB 

TLKO 

TLO. 

1 

TMAGB 

TMAGO 

TMB 

TMO 

TORM 1 

torm 2 

TRB 1 

trb 2 

Tsep 

Tswing 


torques acting at end points of link element 

torque vector on booster c. g. produced by ram mechanism 1 

torque vector on booster c. g. produced by ram mechanism 2 

total thrust torque vector on booster c. g. 
total thrust torque vector on orbiter c. g. 
total thrust vector on booster 
total thrust vector on orbiter 

torque vector at booster ith link attach point from link 
rotations 

total torque on booster c. g. from link forces 

total torque on orbiter c. g. from link forces 

torque vector at orbiter ith link attach point from link 
rotations 

magnitude of booster thrust 
magnitude of orbiter thrust 
total torque vector about booster c.g. 
total torque vector about orbiter c. g. 

torque vector on orbiter c. g. produced by ram mechanism 1 

torque vector on orbiter c. g. produced by ram mechanism 2 

torque at ram 1 attach point on booster 

torque at ram 2 attach point on booster 

release time of separation sequence 

time to initiate swing phase of separation sequence 

time 

displacements at end points of link element 


U1,U2, U4, U6 



U7,U8, U10, U12 
VEB. 

J 

VEO. 

J 

XI, YI, ZI 
A0BL 

C 

0beT , 0OEr 

J 3 
0BL. 

l 

0OL. 

X 

4 B 
40 

[4>B],[<t>0] 

<pBE., <pOE. 

J J 

^B, 0B, (pB 

VB ,0B„, (pB n 
c c c 

^BE., 0BE. 

J J 

^L., <t>L lt 0L. 

^ O, 00, (pO 

^0E„ 0OE. 

3 3 

cdB 


rotations at end points of link element 
vector from booster c. g. to ith booster engine 

vector from orbiter c.g. to ith booster engine 

error signal integrator gains 

link- swing angle 

percent of critical damping 

pitch angles of booster and orbiter engines with roll control 
included 

Euler angle from booster to ith booster link coordinate 
system 

Euler angle from orbiter to ith orbiter link coordinate system 

constant used in ram force relation 
generalized displacements of booster 
generalized displacements of orbiter 
booster and orbiter modal matrices 

hypothetical roll angles of booster and orbiter jth engines 

Euler angles of inertial to booster coordinate system 
booster commanded Euler angles 

gimbal angles of jth booster engine 

relative rotation of end points of ith link member 

Euler angles of inertial to orbiter coordinate system 
gimbal angles of jth booster engine 

angular velocity vector of booster 
angular velocity of orbiter 
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Modifying Symbols: 


B 

refers to 

B., BE. 
3 3 

refers to 

BL. 

refers to 

1 


I 

refers to 

i 

refers to 

3 

refers to 

O 

refers to 

OE. 

refers to 

3 


OL. 

refers to 

l 


X, Y,Z 

body axes 


booster quantity 
jth booster engine 

ith link relative to booster 

the initial values of various quantities 

ith link member 

engines on each vehicle 

orbiter quantity 

jth orbiter engine 

ith link relative to orbiter 

i of either booster or orbiter 


Operators: 

IT 1 

[J T 

0 

n 


inverse of a matrix 


transpose of a matrix 

j / \ 

first derivative with respect to time, 


second derivative with respect to time, 



GENERAL DESCRIPTION OF THE STAGING 
MECHANISM AND THE STAGING SEQUENCE 


The Staging Mechanism 

The link-ram separation system, which is modeled in this thesis, is divided into 
two identical sections. The relative location of these sections is illustrated in fig- 
ure 1(a), and a perspective view of one section is shown in figure 1(b). Each section is 
composed of two standoff links, two parallel swing links, and a ram mechanism 
(fig. 1(b)). The standoff links in both sections and the swing links in the forward sec- 
tion carry all loads between the two vehicles prior to the staging maneuver. The 
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swing-link members in both sections provide a constraint between the booster and 
orbiter during the swing phase of the separation maneuver, and the ram mechanisms 
produce a separation force between the vehicles during this phase. The standoff and 
swing links are pinned parallel to the Y-axis of each vehicle. The ram mechanisms 
are pinned parallel to the Y-axis of the booster and have ball-joint attachments on the 
orbiter. 


The Staging Sequence 

The nominal staging maneuver for the shuttle configuration being modeled (ref. 1) 
occurs at an approximate altitude of 68 885 meters, and the staging velocity is approxi- 
mately 2743 m/sec. The staging sequence is initiated by throttling down the booster 
engines and simultaneously throttling up the orbiter engines as illustrated in figure 2. 
When both vehicles are at the desired thrust levels, the swing phase of the staging 
sequence begins. At this time, the forward and aft standoff links (fig. 1(b)) are pyro- 
technically released from the orbiter, and the ram mechanisms are simultaneously 
activated. The standoff links are retracted into the booster for the duration of the 
flight. During the swing phase of the staging sequence, the orbiter moves forward rel- 
ative to the booster. This relative motion occurs primarily because the thrust-to- 
weight ratio of the orbiter is greater than that of the booster, but the ram forces also 
contribute significantly to the relative velocity between vehicles. 

The duration of the swing phase is scheduled according to time or to link- swing 
angles, whichever achieves its preset maximum value first. At the end of the swing 
phase, the swing links are pyrotechnically disconnected, and the ram mechanisms are 
mechanically disconnected frcm the orbiter. Both the swing links and ram mechanisms 
are drawn back into the booster for the duration of the flight. 


GENERAL DISCUSSION OF THE MATHEMATICAL MODEL 


The formulation of the mathematical model for computing the booster/orbiter 
staging dynamics starts with the determining of the governing equations of motion of 
each vehicle. In this model, each vehicle is given the full six degree of freedom, that 
is, three degrees in translation and three degrees in rotation. The governing transla- 
tional equations of motion are obtained from Newton’s second law, and the governing 
rotational equations are from Euler's moment equations. 

The total forces and torques on each vehicle that are used in the governing equa- 
tions were assumed to be affected by booster and orbiter thrust, interaction link dy- 
namics, ram- mechanism dynamics, and structural dynamics. To facilitate calculation 
of the various forces acting on the booster and orbiter, multiple coordinate systems are 
used in this mathematical model. These coordinate systems define the orientation of 
the booster, booster engines, orbiter, orbiter engines, and link elements. 

To calculate the link forces, each link member is idealized as a beam element. 
The stiffness characteristics for the links are developed, and the forces from the link 
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dynamics are calculated from the linear and angular deflections of the ends of each link. 
The link forces are coupled with the structural response of each vehicle by the vehicle 
response at the link attach points. 

The structural response of each vehicle is calculated from the respective modal 
characteristics of the vehicle and from the thrust, link, and ram forces. This response 
is, in turn, used to calculate vehicle displacements at link, ram, and engine attach 
points. 

The ram mechanism forces are governed by the magnitude of the force generated 
by the ram and by the orientation of the ram between the booster and orbiter. 

The engine thrust vector on each vehicle is a function of the thrust magnitude and 
the thrust vector control (TVC) logic. The TVC logic determines the engine gimbal 
angles based on desired attitude, and the thrust magnitude is controlled by desired 
thrust scheduling. 

This mathematical model was incorporated into a computer program so that the 
staging sequence could be simulated as a function of time. All differential equations 
were numerically integrated in the computer program using an Adams-Moulton pre- 
dictor corrector with a Runge-Kutta starter. The general sequencing of the mathemat- 
ical model is illustrated by the computer program flow chart given in figure 3. The 
formulation of the mathematical model of the staging sequence is discussed in more 
detail in the following chapters. 


COORDINATE SYSTEMS AND TRANSFORMATIONS 


The various coordinate systems used in the mathematical model are illustrated 
in figures 4 and 5. All coordinate systems used are orthogonal, right-hand systems. 
The booster and orbiter systems are coincident with the inertial system at the beginning 
of the swing phase of the staging sequence. The engine and link systems on each vehi- 
cle are oriented relative to that vehicle. 


Inertial Coordinate System 

The inertial coordinate system is fixed with respect to booster and orbiter motion. 
The system is used to define the position and orientation of each vehicle by the position 

vectors SB and SO and the Euler angles i^B, 0B, (pB, and iK), 90, (pO. 

The order of Euler angle rotation from the inertial to booster and orbiter coor- 
dinate systems is about the Z-, Y-, X-axes of the vehicle (fig. 6), and the coordinate 
transformation matrix that corresponds to these angles is defined as follows. 


i 



(i 

T 

I 

GIB(I,J) 

It 

K 


m — 

(k 


( 1 ) 
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where i, j, and k are unit vectors along the respective system axes, and the ele- 
ments of [GIB] are 

GXB( 1 , 1) = COS(eB) • COS(i//B) 

GIB(1, 2) = COS(eB) • SIN(i^B) 

GIB( 1 , 3) = -SIN(eB) 

GIB(2, 1) = -COS(0B) • SIN(^B) + SIN(4>B) • SIN(eB) • COS(i//B) 

GIB(2, 2) = COS(0B) • COS(^B) + SIN(</>B) • SIN(eB) • SIN(i//B) (2) 

GIB(2, 3) = SIN(<£B) • COS(eB) 

GIB(3, 1) = SIN(<£B) • SIN(i//B) + COS(0B) • SIN(eB) • COS(<//B) 

GIB(3, 2) = -SIN(0B) • COS(^B) + COS(<£B) • SIN(eB) • SIN(i//B) 

GIB(3, 3) = COS(tf>B) • COS(eB) 


Also, 


i 




i 

I 


GIO(I, J) 


J 

k, 

lo = 

- 


k 


(3) 


where the elements of [GIO] are defined in the same manner as [GIB] using the orbiter 
Euler angles. Because these Euler angle transformation matrices are orthogonal, then 


[GIB] -1 = [GIB] T 


from which 


T 


I i 



| 

T j 

I j / 

= 

GIB(I,J) 


T 

(kj 

i. 


1 

k) 


(4) 


(5) 
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The Euler angle rate transformation matrix is used to relate vehicle angular ve- 
locities to the Euler angle rates, that is, 


\jsB 



Ubx 

eB 

_ 

GBEA(I, J) 

Lby 

0B 



jffiBZ j 


where the elements of [GBEA] are defined as follows. 

GBEA(1, 1) = 0 

GBEA(1, 2) = SIN(0B)/COS(eB) 

GBEA(1, 3) = COS(0B)/COS(eB) 

GBEA(2, 1) = 0 

GBEA(2,2) = COS(0B) (7) 

GBEA(2, 3) = -SIN(0B) 

GBEA(3, 1) = 1 

GBEA(3, 2) = TAN(eB) • SIN(^B) 

GBEA(3, 3) - TAN(eB) • COS(<£B) 


Likewise, the orbiter Euler angle rates are related to the orbiter angular velocities by 


\p o 




ooOXj 

eO 

= 

GOEA(I, J) 


o)OY 

00 




|cdOZ j 


where the elements of [GOEA] are defined in the same manner as [GBEA] using the 
orbiter Euler angles. 


Engine Coordinate Systems 

The engines on both booster and orbiter are positioned at the engine gimbal point 
and are oriented relative to the respective vehicle coordinate system. Each engine is 
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restricted to rotation about two axes because rotation about the roll (X) axis would not 
produce vehicle control. Thus, the coordinate transformation matrix from the booster 
to the jth booster engine is 


i \ 


| 

i J 

X / 

GBE(I,J) 

' 

T I 

J j 

k 1 



It) 

' ' BE. 





(9) 


where the elements of [GBE] . are as follows. 

J 

GBE(1,1). = COS(eBE.) • COS(<//BE.) 
GBE(1,2). = COS fe BE 
GBE(1,3). = -SIN (e BE.) 

GBE(2, 1). = -SIN^/BE.^ 

GBE(2,2). = COS^BE.) 

GBE(2, 3). = 0 

GBE(3,1). = SIN (e BE.) • COS^/BE.) 
GBE(3,2). = SIN (e BE.) • SIN^BE.) 
GBE(3,3). = COS(eBE.) 



Also, 




IT) 

GOE(I, J) 

It | 


(k) 

1 


( 10 ) 


(ID 


where the elements of [GOE] . are defined in the same manner as [GBE] . using the 
orbiter engine Euler angles. 
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Link Coordinate Systems 


The link coordinate systems are located at the booster and orbiter attach point of 
each link. The system at each end is oriented with respect to the coordinate system 
of the vehicle to which it is attached. Because each link member is pinned parallel to 
the Y-axis of the attached vehicle, the link coordinate system is free to rotate only 
about the Y-axis of that vehicle. Thus, the coordinate transformation matrix from 
booster to the ith booster link coordinate system is 


m 



i 

It! 

GBL(I, J) 


I 

I k /BL. 

- 

1 

[k 


where the elements of [GBL] . are 

GBL(1, 1). = COS(eBL.) 
GBL(1,2). = 0 
GBL(1, 3). - -SIN(eBL.) 
GBL(2, 1). = 0 
GBL(2,2). - 1 
GBL(2,3). = 0 
GBL(3, 1). = SIN(eBL.) 
GBL(3, 2). = 0 
GBL(3,3). = COS(eBL.) 

Also, 


i 




i J 

I 

k 

GOL(I, J) 


J | 

k 

i J 

' OL. 



! k lo 


( 12 ) 


(13) 


(14) 


where the elements of [GOL] . are defined in the same manner as [GBL]. using the 
orbiter link rotation qOL^. 1 1 


14 



VEHICLE EQUATIONS OF MOTION 


The equations of motion for the booster and orbiter are obtained from Newton’s 
second law for translation and from Euler’s moment equations for rotation. 

The elementary equations of translation for each vehicle in component form are 


(SBX) 1 

FSMBX] 

l .. 1 l ' 

| f 

SBY ■ MBT ; 

FSMBY 

(SBZ) 1 

[fsmbz) 
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SOZ 

1 1 

[FSMOZ 


(15) 


(16) 


where FSMBX, FSMBY, FSMBZ and FSMOX, FSMOY, FSMOZ are the compo- 
nents of the total external forces on each vehicle, respectively. 

Euler's moment equations for rotational motion for the booster (refs. 2 and 3) are 
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(17) 


where TMBX, TMBY, TMBZ are the components of the total torque vector on the 
booster c. g. 
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The elements of [INRB] are 


INRB(1, 1) = EXXB 
INRB(1,2) = -IXYB 
INRB(1, 3) = -EXZB 
INRB(2, 1) = -IXYB 

INRB(2,2) = IYYB (18) 

INRB(2, 3) = -IYZB 

INRB(3, 1) = -IXZB 

INRB(3, 2) = -IYZB 

INRB(3, 3) = IZZB 


Similarly, Euler’s moment equations for the crbiter are 
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(19) 


where the elements of [INRO] are defined in the same manner as [INRB] using the 
orbiter inertias. 

The translational equations of motion can be integrated with respect to time to 
yield the booster and orbiter inertial velocities. Likewise, these velocities can be inte- 
grated to yield the corresponding displacements of each vehicle. 
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and 


also 


and 
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The rotational equations (17) and (19) can be solved for <dB and mO, respectively. 
These vehicle angular accelerations can be integrated with respect to time to yield 
vehicle angular velocities, that is, 
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(24) 


(25) 
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The vehicle angular velocities can be transferred into Euler angle rates (ref. 2) by the 
transformation matrix of equations (6) and (8), that is, 
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(26) 


(2.7) 


The Euler angle rates can then be integrated with respect to time to yield the booster 
and orbiter Euler angles, that is, 
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and 
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(29) 


Thus, the position and orientation of the booster and orbiter at any time can be deter- 
mined by knowledge of the external forces acting on the respective vehicles. 

As previously mentioned, these external forces are affected by link and ram dy- 
namics and the thrust of each vehicle. The total force vectors on each vehicle are 


FSMB = FLKB + BRMV + THRVB 


(30) 
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and 


FSM$ = FLKO + ORMV + THRVd (31) 

where FLKB and FLKO are total link forces, BRMV and ORMV are total ram 
forces, and THRVB and THRVO are total thrust forces on the booster and orbiter, 
respectively. 

The corresponding total torques acting on each vehicle are 


TMB = TLKB + TBRM + THMB (32) 

and 

TMO = TLKO + TORM + THMO (33) 

where TLKB and TLKO are total link torques, TBRM and TORM are total ram 

torques, and THMB and THMO are total thrust torques on the booster and orbiter 
centers of gravity, respectively. 

The development of the various forces and torques acting on the booster and 
orbiter is discussed in the following sections. 


LINK FORCES 


The forces generated in each link element are calculated from the translational 
and rotational displacements at each end of the link element and from the corresponding 
stiffness characteristics of that element. 

The translational displacements of each link element are calculated as follows. 
The vector SEPV. (fig- 7) between the booster ith link attach point and the correspond- 
ing orbiter ith attach point is calculated from vector loop equations. A vector from the 
booster c. g. to the orbiter ith link attach point can be written in two forms. Because 
the two equations define the same point from the same origin, they can be equated, 
that is, 


POSB. + SBD. + SEPV. = SBO + POSO. + SOD. 

ill ii 


(34) 
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where PQSR and POSO. define the ith link attach point on each vehicle relative to its 

e.g. The vectors SBD. and SOD. are the structural displacements at the ith link 

attach point on each vehicle, respectively. The vector SBO defines the orbiter c. g. 
relative to the booster c. g. and is calculated from the inertial position vector of each 
vehicle, that is, 


SB + SBO = SO 


or 


SBO = SO - SB 


(35) 


(36) 


Equation (36) is substituted into equation (34), and solving for SEPV. yields equa- 
tion (37). 1 

SEPV. = SO - SB + POSO. + SOD. - POSB. - SBD. (37) 

l 1111 ' ' 


that defines the position of the ith link attach point on the orbiter with respect to the 
booster ith link attach point. This vector can be expressed in the link coordinate sys- 
tem, and then the deflections of that link can be directly calculated from the vector 
components, that is, 


DLX. = SEPVX. (t) - SPXI. 

l l w l 


(38) 


and 


DLY. = SEPVY. (t) 

l l ’ 


(39) 


The link deflection in the ZBL and ZOL 
link is pinned about the YBL and YOL axis. 


direction (fig. 5) is zero because the 
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/ 


The rotations of the end points of each link element are calculated as follows. 
Vectors in the booster, booster link, orbiter, orbiter link, and inertial coordinate sys- 
tems are related according to equations (1), (12), (13), and (14), that is, 
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(41) 


(42) 


(43) 


Equations (40) and (41) can be written as 
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and 
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Equation (45) can be substituted into equation (44) to yield the coordinate trans- 
formation from the booster link to the inertial system, that is, 
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BL. 


This relation (eq. (46)) can then be substituted into equation (42) to yield equa- 
tion (47) as given below. 
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Equation (47) can then be substituted into (43) to yield the coordinate transforma- 
tion from booster link to orbiter link coordinate systems, that is, 
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where 


[GBOL(I,J)] i = [GOL(I,J)].[GIO(I,J)][GIB(I,J)] T [GBL(I,J)] 1 T 


(50) 


Now the rotations of ith link at the orbiter attach point (link end II) with respect to its 
booster attach point (link end I) can be solved for by the ratio of certain elements of 
[GBOLj.. The elements of [GBOL]. for Euler angle rotation order Z, Y, X are 

given below. 


GBOL(l, 1). = COS(eL.) • COS(^L.) 

GBOL(l,2). = COS(eL.) * SIN^L.) 

GBOL(l,3). = -SIN(eL.) 

GBOL(2, 1). = -COS^L.) • SIN^L.) + SIN SIN (e L.)COS^L.) 
GBOL(2, 2). = COS(0L.)COS(i//L.) + SIN^L^SIN^e LjSIN^L.) 
GBOL(2, 3). = SIN^L.) • COS(e L.) 

GBOL(3, 1). = SIN^L^SIN^L.) + COSp>L.)SIN(e L.)COS^L.) 
GBOL(3,2). = -SIN^L^COS^L.) + COSp>LjSm(e L.JSIN^L.) 
GBOL(3, 3). = COS(4>L.)COS(eL.) 


(51) 


from which 

0L. = TAN' 1 (GBOL(2 , 3).) /(GBOL(3, 3).) 
j//L. = TAN" 1 (GBOL(l,2).y(GBOL(l,l).) 


The link stiffness characteristics are calculated by idealizing each link member 
as a beam element. The boundary conditions for each link can then be applied to the 
general stiffness matrix of a beam element as given in reference 4 to yield the stiffness 
matrix of each link (appendix A). The translational and rotational displacements of the 
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ith link end II relative to end I can be multiplied by the corresponding link stiffness 
matrix to yield forces at the orbiter attach point, that is, 


(FLOX.'l 


"AE/£ 

0 

0 

0 


ft) LX.) 

l 

FLOY. 


0 

12EI/1 3 

o 

-6EI/I 2 


l 

DLY. 

\TLOX. 

1 

> B 

0 

0 

jG/e 

O 

< 

(PL. > 

TLOZ. 

1 ij 


0 

-6EI/1 2 

0 

4EI/I 


r i j 


V v L 1 ) 


and the reaction forces on the booster are given in equation (54). 
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( 53 ) 


(54) 


The total link forces and torques acting on each vehicle are the vector sum of the in- 
dividual link forces and torques, that is, 
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RAM MECHANISM FORCES 


(55) 


(56) 


The total ram mechanism interaction forces between the booster and orbiter are 
a function of the forces produced by activation of the mechanisms and the forces be- 
cause of structural deflection of each mechanism. 
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The force magnitude produced by the activation of each ram mechanism was as- 
sumed to have a predetermined time history as illustrated in figure 8. The force equa- 
tion of each ram is 


FRAM = FRM(l - e“ Xt ) ( 57 ) 


where FRM is the maximum ram force and X is a positive constant. The force vector 
from the above ram force magnitude is calculated by multiplying this magnitude times 
a unit vector along the line of action of the ram. The orientation of each ram is calcu- 
lated in the same manner as the link members. The loop equations that define the ram 
attach point on the orbiter relative to the booster c.g. can be equated. Therefore, for 
ram mechanism 1 (forward mechanism) 


PBSRj + SBD r1 + SRMj = SBO + POSR 1 + SOD R1 (58) 

Equation (58) can be solved for the vector SRMj, which defines the line of action of the 
ram force. Thus, 


SRMj = SBO + POSRj + SOD R1 - PBSRj - SBD Rl 


(59) 


The forces on the booster and orbiter because of activation of ram 1 are 


BRMV 1 = -FRAM * SRM^ / 


SRM, 


(60) 


and 


ORMV 1 = FRAM • SRM-JJsRM^ 


(61) 


The corresponding torques on each vehicle are 

'l 1 =(PBSR 1 + SBD r Jx BRMV ] 


TBRM, 


(62) 
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and 


TORMj = ^POSRj + SOD R1 )x ORMY 1 


(63) 


The structural forces produced by the ram mechanisms are a function of the ram 
stiffness characteristics and the structural displacements of the mechanisms. As pre- 
viously mentioned, the ends of the ram mechanisms are pinned parallel to the Y-axis 
of the booster and attached as a ball joint on the orbiter. The lateral (YB-direction) 
displacements of the orbiter ram attach point relative to its booster attach point are 
calculated from equation (59), that is, 


DLY r1 = SRMYj - SYIj (64) 


where SYI^ is the unstroked length of ram 1. 

This deflection is multiplied by the lateral stiffness of the ram mechanism to 
yield the lateral force at the orbiter attach point. 


FRY. = -KRL, • DLY_. 
1 1 R1 


(65) 


The reaction force at the booster attach point is simply 


FRBYj = -FRYj 


( 66 ) 


The corresponding torque at the booster attach point because of the lateral deflection 
force is 


TRBj 


SRM 1 


• FRYj 


(67) 


Because of ram deflection, the lateral force acts along the YB-axis, and the torque acts 
about the XB-axis (fig. 1). Thus, the force and torque vectors caused by ram mecha- 
nism deflections are 


FRBj = FRBYj j R 


( 68 ) 
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and 


TRBj = TRBj ig 


(69) 


Because of ram deflections, the force vector on the orbiter is calculated by transfer- 
ring FRY^ from the booster to orbiter coordinate system. Thus, from equations (1) 

and (3), the force on the orbiter because of ram deflection is given by equation (70) 
below. 
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(70) 


The forces from ram mechanism 2 (rear mechanism) are calculated in the same man- 
ner as for ram mechanism 1. 

The total ram mechanism forces acting on the booster and orbiter are 


BRMV = BRMVj + BRMVg + FRB^ + FRB, 


(71) 


and 


ORMV = ORMVj + ORMV 2 + FROj + FROg 


(72) 


The corresponding torques are 

TBRM = TBRMj + TBRM 2 + TRBj + TRB 2 +(pBSR 1 + SBD R1 ) x FRBj 

+ (PBSR 2 + S^ R2 ) X FRB 2 (73) 

and 

TORM = TORMj + TORMg + (pOSRj + SOD R1 ) x FR^ + (pOSR 2 + X ™2 

(74) 
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THRUST FORCES 


The thrust vector from each engine on the booster and orbiter was calculated 
from the thrust magnitude and from the engine gimbal angles. The thrust magnitude 
on both the booster and orbiter was scheduled to achieve desired thrust- to- weight ra- 
tios. The thrust scheduling of each vehicle is given in figure 2. The engine gimbal 
angles are calculated from the TVC logic, which is presented for planar rotations in 
reference 5 and is extended to cover full vehicle rotations. 

The TVC logic is described for the booster engines as follows. The booster atti- 
tude deviations from the commanded (desired) attitude about each axis are 


DTBl = - »^BC 

DTB2 = 0 B - 9 BC 
DTB3 = (pB - 4>BC 


(75) 


These attitude deviations are transferred into the booster coordinate system so that 
they can be combined with the booster angular velocity components to calculate the 
booster error signals. The attitude deviations are calculated by differencing Euler 
angles. These deviations are small; thus, they are assumed to be vector quantities. 
The deviations are transferred into the booster system using the inverse of equation (6). 
The attitude deviations expressed in the booster coordinate system are 
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(76) 


The error signals of the booster are calculated as a weighted sum of the attitude devia- 
tion and the angular velocity components. The weighting factors are called gains, and 
they reflect the desired weighting of the component by which they are multiplied. 
Therefore, 


EBX = ETBX + RX • 

03 BX 


EBY = ETBY + RY • 

«3 BY 

(77) 

EBZ = ETBZ + RZ * 

03 BZ 



28 



These error signals are integrated to yield the error integrals, that is, 

EBXI f EBX EBXIj 

EBYI = EBY dt + EBYIj (78) 

EBZI I EBZ EBZL 

i J \ / i l / 

The gimbal angles and rates of each booster engine required to correct the booster 
attitude deviations are calculated from the error signals and error integrals. The gim- 
bal angle commands for the jth booster engine are calculated as a weighted sum of the 
error signals and the error integrals, that is, 

GBCX. = DX Bj (EBX + XI Bj • EBXl) 

GBCY. = DY Bj (EBY + Ylg. * EBYl) (79) 

GBCZ. = DZg^EBZ + ZI Bj • EBZlj 

Because the engine cannot move instantaneously to this commanded angle, the gimbal 
rate command is calculated. These commands for the jth engine are 

GBXD. = CXg^GBCX. - 0BE.j 

GBYDj = CY b .(GBCY. - 0 BE. j (80) 

GBZD. = CZg.^GBCZ. - ^BE.j 

The gimbal rates are limited by the physical size of the engine and the mechanism that 
gimbals the engine. The gimbal angles of the jth engine can be calculated from the 
gimbal rates and the initial gimbal angles, that is, 

0BE\ f GBXD] 4 >BEj 

, q BE > = GBYD dt + 0 BEj 

il/BE . GBZD . 1//BE. 

' j J \ 1 ' L 
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The engine roll angle <£BEj can be combined with engine pitch or yaw angles to pro- 
duce roll control. For the case of two side-by-side engines, the roll angles can be 
combined with the pitch angles for the roll control, that is, 


0 BE ' = e BEj + KRRj • (p BEj 
0 BEg = 0BE 2 - KRR 2 • (pBE 2 


(82) 


These equations represent one engine pitching up while the other pitches down to create 
the desired roll torque. 

Thus, the gimbal angles of each engine are known, and, assuming the engine 
thrust acts along the centerline of the engine, the engine thrust vector can be calcu- 
lated, that is, 
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The engine gimbal angles of the orbiter are calculated in the same manner as 
those of the booster; thus, the orbiter thrust vector is 
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The total booster and orbiter thrust are 
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The corresponding c.g. torques are 


THMB 


Ml . 
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(87) 


and 
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1V14 / ^ ^ V 

THMO = ( VEO j X THRVO j) 


( 88 ) 


STRUCTURAL DYNAMICS 


The structural response of each vehicle can be calculated from the external load- 
ing and the modal characteristics of that vehicle. The external loads on each vehicle 
are the link forces, ram forces, and thrust forces. The modal charactenstics of ea« ch 
vehicle are calculated from their respective structural models. The modes and fre- 
quencies of a multi degree-of-freedom model of both booster and orbiter can be devel- 
oped as described in reference 6. 

For this mathematical model, the structural response of the vehicles is included 
to examine its effect on the interaction link forces between vehicles. Because the 
forces generated in each link member are a function of the position and orientation of 
SSS, the structural response of each vehicle at the link attach points is 

of interest. 

The structural displacement at the link attach points can be related to the struc- 
tural response of each vehicle by the following equations. 


{SBD} = [®B] |4B[ 

(89) 

{SOD} = [®0] {40} 

(90) 


where [0B] and [00] are the 
£B and 40 are the generalized 
displacements at specific points 


"free-free" modal matrices of the booster and orbiter, 
displacements, and SBD and SOD are the structural 
of interest. 
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The generalized displacements are calculated in the following manner. The gen- 
eralized acceleration of the booster as developed in reference 6 can be written as 


l*Bt ={MBG] [$B] T jFBEf -2 c[bMg]|/b} - |^BMG^ j£B} (91) 

where [MBG] is the generalized mass matrix and [BMG] are the natural frequencies 
for the booster structural model. The constant £ represents the percent of critical 
damping, and |FBE} is the vector of external forces acting on the booster. The gen- 
eralized acceleration of the orbiter is of the same form as equation (91). 

The generalized accelerations of each vehicle can be integrated to yield the cor- 
responding generalized velocities and displacements, that is, 


\£ b \ = /||Bldt + ||Bj| 

\£o\ = J{£)ldt + (92) 


and 


|£B} = /|4B}dt + {iBj 

= /{&>fdt + ^OjJ. (93) 


The generalized displacements of equation (93) can be used in equations (89) and (90) to 
yield structural displacements of each vehicle at the desired node points (fig. 9). 


COMPUTATIONS 


As noted in the introduction, the mathematical model is incorporated into a com- 
puter program so that the staging maneuver can be simulated. The computer program 
is described in appendix B. The primary input parameters used in the staging simula- 
tions are divided into two categories: those that were fixed for all simulations and 
those that were varied between simulations. 


Fixed Input Data 

The masses, inertias, and c. g. locations of each vehicle are given in table I. 
These quantities were assumed to be fixed with respect to each vehicle over the dura- 
tion of the separation maneuver. The inertial coordinate system was positioned at the 
initial booster c. g. location, and the c. g. of the orbiter is defined relative to this 
initial position. 
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The position vectors of the link and ram mechanism attach points relative to the 
c. g. of each vehicle are given in table H. Also, the position vectors of the booster and 
and orbiter engines relative to their respective centers of gravity are presented in 
table II. For the simulations presented in this report, the booster and orbiter were 
each assumed to have one engine with equivalent thrust of the multiple engines of the 
respective vehicles described in reference 1. The TVC gains for each vehicle are 
given in table in. The booster engine TVC logic was used during the swing phase, but 
the orbiter engine gimbal angles were constant over this period. 

The vehicle structural data used in the simulations were very limited. Only lim- 
ited data were available because of the preliminary stage of development of the Space 
Shuttle Program. The numbering of the orbiter node points is illustrated in figure 9, 
and the corresponding modal matrix for the orbiter is given in table IV. The orbiter 
natural frequencies and generalized masses also are given in table IV. No booster 
structural data were used because no data were known to exist for the configuration 
being simulated. 


Variable Input Data 

Certain quantities were varied between simulations so that the effects of these 
variations could be examined. These quantities include vehicle thrust scheduling, ram 
forcing function, and link elastic properties. 

Case 1 of the simulations represented the nominal case for the link-ram system. 
The variable input data, coupled with the fixed input data, represent the baseline input 
for the nominal separation sequence. The subsequent cases are simulated with pertur- 
bations in the variable input data of case 1. The thrust scheduling for case 1 is pre- 
sented in figure 2, and the ram forcing function for this case is presented in figure 8. 
The link stiffness characteristics for this system are given in table V. No elastic 
modes were used in cases 1 to 4. 

Case 2 was simulated with a 100-percent increase in the stiffness of each link 
member. 

Case 3 was simulated with different ram mechanism force characteristics. The 
ram force was simulated as a step function (X = °°) and for a lower than nominal onset 
rate (X = 8) as shown in figure 8. 

Case 4 was simulated with the orbiter at full thrust at the beginning of the swing 
phase, as opposed to approximately 67 percent of full thrust in the nominal case. 

Case 5 represented the nominal separation sequence with the first orbiter mode 
included and the nominal separation sequence with the first 10 orbiter modes included. 
These modes correspond to the lowest natural frequency of the orbiter. 
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RESULTS 


Results are presented for the five cases of the separation maneuver that were 
simulated. The parameters that are presented include swing-link loads, swing-link 
rotation angle, booster and orbiter angular velocities, relative c. g. translational ve- 
locities, and booster engine gimbal angle. 

The load time histories in the forward and rear swing links for case 1 are pre- 
sented in figures 10 and 11. The loads in the forward swing links (fig. 10) oscillate 
between compression and tension during the swing phase of the staging sequence. The 
loads in the rear swing links (fig. 11) remain primarily in tension during the swing 
phase. Tension loads in both the forward and aft links produce negative torques on the 
booster c. g. because the attach points of these links are forward of the booster c. g. 
(table n). Tension loads in the forward and aft links also produce negative torques on 
the orbiter because of the position and orientation of the links relative to the orbiter 
c. g. (table II). 

The angular velocities about the Y-axis of each vehicle for case 1 are presented 
in figures 12 and 13. The X- and Z -components of angular velocity were zero in all 
cases simulated. In all cases, both the booster and orbiter build up a negative (nose 
down) angular velocity during the swing phase. The booster engine gimbal angle is 
presented in figure 14. This gimbal angle is decreasing throughout the swing phase 
because the booster engine is trying to correct the negative booster angular velocity. 

The swing angles of both the forward and rear swing links are approximately equal. 

The forward link- swing angle (fig. 15) is plotted relative to its initial angle, which is 
17.7°. The engine gimbal angle and the link- swing angle are shown only for case 1 
because only small changes in these values existed in the subsequent cases. 

The relative axial (XI- direction) and normal (Zl-direction) velocities between the 
booster and orbiter centers of gravity are presented in figures 16 and 17, respectively. 
These linear velocities, coupled with the angular velocities, can be used to calculate the 
relative velocity between any points on the booster and orbiter. 

The effects of increased link stiffness on the staging dynamics were examined in 
case 2. The forward and rear swing-link loads are presented in figures 18 and 19. The 
frequency of oscillation in both the forward and rear swing-link loads is about 25 percent 
higher than nominal. The higher frequencies would be expected for the increased link 
stiffnesses. The maximum value of the forward link loads (fig. 18) are about the same 
as nominal, but the peak tension loads in the rear swing links (fig. 19) are about 15 per- 
cent higher than nominal. The angular velocities for this case are presented in fig- 
ures 20 and 21. The booster angular velocity (fig. 20) is about the same as the nominal 
case, but the orbiter angular velocity (fig. 21) is more oscillatory, and its peak value 
is about 10 percent lower than nominal. The relative translational velocities of this 
case are presented in figures 22 and 23. These values are virtually unchanged from 
the nominal case. Thus, the increased link stiffness does significantly change the rel- 
ative motion between the booster and orbiter during the separation maneuver. 

Case 3 was simulated to examine the effects of different ram mechanism forces 
on the separation performance. The step input ram force (A. = °°) represented about a 
20-percent increase in impulse; the lower onset ram force (A = 8) represented a 



30-percent decrease in impulse. The loads in the forward and rear swing links for this 
case are presented in figures 24 and 25. The forward link loads (fig. 24) for the step 
ram force are approximately 35 percent higher than nominal in compression; the rear 
link loads (fig. 25) for this ram force were about 35 percent higher than nominal in 
tension. The maximum loads in both the forward and rear links for the lower onset 
ram force are about the same as nominal. The booster and orbiter angular velocities 
for this case (figs. 26 and 27) build up negative values faster for the step ram force 
than for the lower onset ram force. However, the booster angular velocities (fig. 26) 
at the time of release (Tsep) are approximately equal for both ram forces and are 
approximately 5 percent below nominal. Also, the orbiter angular velocities (fig. 27) 
at release are approximately equal for both ram forces and are approximately 10 per- 
cent below nominal. The relative axial velocity (fig. 28) for the step ram force is about 
15 percent higher than nominal and about 10 percent below nominal for the lower onset 
ram force at release. Similarly, the relative normal velocity (fig. 29) for the step 
ram force is about 10 percent higher than nominal, and it is 10 percent below nominal 
for the lower onset ram force. Thus, the relative motion of the booster and orbiter is 
proportional to the total impulse of the ram-mechanism forces. However, this propor- 
tionality is not one to one because the thrust forces also contribute to the total impulse 
on each vehicle. 

Case 4 was run with higher than nominal orbiter thrust at the beginning of the 
swing phase (Tswing). This higher orbiter thrust level represented a 30-percent in- 
crease in impulse to orbiter during the swing phase. The forward link load time his- 
tory for this case (fig. 30) has peak tension and compression loads that are about the 
same as nominal. Likewise, the maximum values of the rear link loads (fig. 31) are 
about equal to nominal. The booster and orbiter angular velocities (figs. 32 and 33) 
are each about 12 percent higher than nominal at release. The relative axial velocity 
for this- case (fig. 34) is approximately 20 percent higher than nominal at release, and 
the relative normal velocity (fig. 35) is approximately 5 percent higher than nominal at 
release. Thus, the increased orbiter thrust in this case represents an increase in im- 
pulse acting axially on the orbiter, and this increase is reflected in the larger relative 
axial velocity. 

The effects of the orbiter modal data on the separation performance are examined 
in case 5. The forward swing-link loads (fig. 36) for one mode and those for 10 modes 
are similar in magnitude but slightly different in frequency. The maximum tension and 
compression values of these loads are about 15 percent below nominal. The rear link 
loads (fig. 37) for one mode and those for 10 modes are also similar in magnitude and 
slightly different in frequency. The magnitude of these loads is about equal to nominal 
although the maximum tension load for 10 modes is about 8 percent higher than nominal. 
The booster angular velocities are presented in figure 38 and are about the same as 
nominal for one mode and for 10 modes. The orbiter angular velocities (fig. 39) are 
about 15 percent below nominal at release for both one mode and 10 modes. The rela- 
tive axial and normal velocities for this case (figs. 40 and 41) are virtually unchanged 
from nominal for both one mode and 10 modes. Thus, including the orbiter structural 
response does not significantly change the relative motion of the booster and orbiter 
during the separation maneuver. 

The linear and angular velocities of the booster and orbiter were used to calculate 
the separation distance between points on each vehicle during the separation maneuver. 
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The separation distance at the various points was used as a check for possible recon- 
tact between the vehicles. In all simulated cases, no recontact between the booster 
and orbiter occurred. 

The stresses resulting from both the tension and compression loads in the swing 
link members were well below the proportional limit stress of these members for all 
cases simulated. 


CONCLUSIONS 


The foregoing discussion of the staging dynamics of a space shuttle configuration 
using a link-ram staging mechanism leads to the following conclusions. 

1. For all cases simulated, the translational and rotational motion of the booster 
and orbiter during the separation maneuver is sufficient to ensure no recontact between 
the vehicles. 

2. The loads generated in the swing-link members during the swing phase of the 
staging sequence produce stresses that are well below the proportional limit stress of 
these members. 

3. Variations in the stiffnesses of the swing-link members change the maximum 
loads in these members but do not significantly change the relative motion of the vehicle 
during the staging maneuver. 

4. Variations in the total impulse from the orbiter thrust and from the ram mech- 
anism forcing function produce appreciable change in the relative motion of the booster 
and orbiter during the staging maneuver. 

5. Including the orbiter structural response in the staging simulations produces 
only a small deviation in the link loads and a negligible change in the relative motion of 
the vehicles. 

The following recommendations are made concerning the further analysis of the 
staging maneuver. 

1 . Booster heating caused by orbiter engine plume impingement should be 
evaluated. 

2. The effects of plume impingement forces on the staging dynamics should be 
examined. 

3. The nonsymmetric loading on each vehicle because of deviations in link re- 
lease times should be examined. 
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4. The effects of fuel slosh on the staging dynamics should be evaluated. 

5. The performance of the link-ram staging system under high aerodynamic 
load conditions should be evaluated. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, October 11, 1972 
986-15-31-06-72 
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TABLE I. - VEHICLE MASS PROPERTIES 








TABLE II. - LINK, RAM, AND ENGINE LOCATION 



Axis 


Axis 


Vector ' 

X 

Y 

Z 

Vector 

X 

Y 

z 



Link attach points, m 



POSB 1 

30. 68 

-1. 777 

-6. 65 

POSOj 

3.03 

-1. 777 

2. 17 

posb 2 

30.68 

-1. 777 

-6. 65 

POS0 2 

3.03 

-1. 777 

2. 17 

POSBg 

26. 11 

-1. 777 

-6. 65 

POSOg 

3. 03 

-1. 777 

2. 17 

posb 4 

26. 11 

-1. 777 

-6. 65 

POS0 4 

3.03 

-1. 777 

2. 17 

POSBg 

10.26 

-1. 777 

-6. 65 

POSOg 

-17. 39 

-1. 777 

2. 17 

POSBg 

10.26 

-1. 777 

-6. 65 

POSOg 

-17. 39 

-1. 777 

2. 17 

POSB ? 

5. 69 

-1. 777 

-6. 65 

poso 7 

-17. 39 

-1. 777 

2. 17 

POSBg 

5. 69 

-1. 777 

-6. 65 

POSOg 

-17. 39 

-1. 777 

2. 17 



Ram mechanism attach points, m 



PBSRj 

30. 68 

0 

-6. 65 

POSRj 

3.03 

0 

2. 17 

pbsr 2 

5.69 

0 

-6. 65 

POSR 2 

-17. 39 

0 

2. 17 

1 — 

Engine gimbal points, 

m 



BEV 

-21.56 

0 

-1. 77 

OEV 

-15.26 

0 

-. 648 
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Vehicle 

uain 

Booster 

Orbiter 

RX 

0.7 

2 

RY 

.7 

2 

RZ 

.7 

2 

DX 

2 

2 

DY 

2 

2 

DZ 

2 

2 

XI 

. 1 

.05 

YI 

. 1 

.05 

ZI 

. 1 

.05 

CX 

15 

15 

CY 

15 

15 

CZ 

15 

15 









TABLE IV. - ORBITER MODAL MATRIX, NATURAL FREQUENCIES, AND GENERALIZED MASSES 


Mode 


Node 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

IX 

0. 00907 

-0. 00340 

0. 06036 

-0.04472 

0. 07836 

-0.13010 

-0.13955 

-0. 81576 

0. 65862 

0.17119 

1Z 

-.70784 

-.97249 

. 00272 

. 00894 

.21479 

-.46816 

. 06518 

. 07037 

.10938 

.01714 

2X 

. 00907 

-.00340 

. 06036 

-.04472 

.07836 

-.13010 

-.13955 

-.81576 

.65862 

.17119 

2Z 

-.70784 

-.97249 

. 00272 

. 00894 

.21479 

-.46816 

. 06518 

. 07037 

.10938 

.01714 

3Z 

. 99642 

-.46076 

.11946 

-.04861 

-.85350 

-.71552 

-.73947 

.19871 

-.13052 

-.02099 

4Z 

. 99642 

-.46076 

.11946 

-. 04861 

-.85350 

-. 71552 

-.73947 

.19871 

-.13052 

-.02099 

1 

5X 

-.14138 

.12760 

-.08165 

-.02682 ! 

. 09469 

-.05436 1 

.20714 

. 34208 

-. 33816 

-.02914 

5Z 

1.3894 

-.86968 

. 33203 

-.05869 

-1.1818 

-1. 0365 

-2.1068 

. 50503 

-.27140 

. 27498 

6X 

. 00907 

-.00340 

. 06036 

-.04472 

.07836 

-.13010 

-.13955 

-.81576 

. 65862 

.17119 

6Z 

-.70784 

-.97249 

. 00272 

. 00894 

.21479 

-.46816 

. 06518 

. 07037 

.10938 

.01714 

7Z 

. 99642 

-.46076 

.11946 

-.04861 

-.85350 

-.71552 

-.73947 

.19871 

-.13052 

-.02099 


Notes: (1) = 2. 8697 

(7) co ? = 10. 5549 

(2) cd 2 = 4.8610 

(8) cog = 11. 7859 

(3) co<j = 6.8383 

(9) co 9 = 11.9221 

(4) cd 4 = 7. 5418 

(10) cd^q = 12. 9829 

(5) cc g = 7.9059 

(11) MOG. = 173227. 3 

(6) cd 6 = 8.0210 

(12) i = 1, 10 


Dimensions: (1) MODAL MATRIX = m/m 

(2) co = Hz 

(3) MOG = kg 




TABLE V. - LINK STIFFNESS CHARACTERISTICS 

E = 1. 1 x 10 11 N/m 2 ; G = 4. 068 x 10 10 N/m 2 ; 

8 2 

proportional limit stress = 8. 55 x 10 N/m 


Link 

1 

2 

3 

4 

5 

6 


8 


Properties 


l , m 

A 2 

A, m 

I, 

J, 

4.829 

0.0110 

0.00018 

0.00036 

4.829 

.0110 

.00018 

.00036 

1.473 

.0039 

.0035 

.0069 

1.473 

.0039 

.0035 

.0069 

4.829 

.0064 

. 000035 

.000071 

4.829 

. 0064 

. 000035 

.000071 

1.473 

.0039 

.0035 

. 00070 

1.473 

.0039 

.0035 

. 00070 





















Figure 3.- Computer program flow chart. 







Booster c.g. Booster engine gimbal point 


Note: Ail Y-axes are parallel 
Figure 4.- Illustration of coordinate systems. 


Orbiter attach point 



Booster attach point 

Figure 5.- Link element coordinate systems. 




(a) The rotation of XY-plane about 
Z-axis. 


(b) The 9 rotation of X'Z' -plane about 
Y'-axis. 



(c) The <p rotation of Y"Z"-plane about X"-axis. 
Figure 6. - Euler angle rotations. 







Figure 8. - Raw mechanism force 
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Figure 9. - Orbiter node points where external forces are applied. 
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Figure 11.- Case 1, force in rear swing links 




Booster engine gimbal ar 




Figure 16. - Case 1, relative inertial axial velocity between booster 
and orbiter centers of gravity. 



Relative normal velocity, m/ 



Figure 17.- Case 1, relative inertial normal velocity between booster 
and orbiter centers of gravity. 
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Figure 23. - Case 2, relative inertial normal velocity between booster 
and orbiter centers of gravity. 











Figure 29. - Case 3, relative inertial normal velocity between booster 
and orbiter centers of gravity. 



Force in forward swing links. 


-556.0 x 10 3 



Figure 30. - Case 4, force in forward swing links. 
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-44.48 x 10 3 



Figure 31.- Case 4, force in rear swing links. 
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Figure 33.- Case 4, orbiter angular velocity. 
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Relative normal velocity, m/sec 



Tswing Time, sec Tsep 

Figure 35. - Case 4, relative inertial normal velocity between booster 
and orbiter centers of gravity . 
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Figure 37. - Case 5, force in rear swing links. 
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Relative axial velocity , m/sec 
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1st mode 
1st 10 modes 



Tswing Time, sec Tsep 


Figure 39. - Case 5, orbiter angular velocity. 



Figure 40. - Case 5, relative inertial axial velocity between booster 
and orbiter centers of gravity. 





APPENDIX A 


LINK STIFFNESS CALCULATIONS 


The link stiffness calculations are made by idealizing each link member as a 
beam element. The general stiffness for a beam element that undergoes a full six- 
degree -of -freedom deflection at each end is derived in reference 4. The link element 
boundary conditions can be applied to the general beam stiffness matrix so that the link 
stiffness matrix can be formulated. These boundary conditions are that each link mem- 
ber is pinned at its end points with the pin parallel to the Y^-axis (fig. A-l). The 

deflections and forces that correspond to the link element boundary conditions, as 
illustrated in figure A-l, are related by the link stiffness matrices; that is, 


(Y\ \ 


AE/£ 








fui \ 

F2 


0 

12El/£ 3 




Symmetric 



U2 

T4 


0 

0 

JG/£ 






U4 

1 T6 


o 

6El/£ 2 

O 

4El/£ 




< 

U6 

|F7 

-AE/£ 

O 

O 

O 

AE/£ 



U7 

F8 


0 

-12El/£ 3 

o 

-6El/£ 2 

0 

12El/£ 3 



U8 

T10 


0 

O 

-JG/£ 

O 

0 

0 JG/£ 



U10 

0ri2j 


0 

6El/£ 2 

O 

2El/£ 

o 

-6El/£ 2 0 

4El/£ 


JJ12y 


(A-l) 


where A is the link element cross-sectional area, I is the cross-sectional moment of 
inertia, J is the polar moment of inertia, E is the modulus of elasticity, G is the 
shear modulus of elasticity, and £ is the link element. 

For simplification, the deflections of link end II (fig. A-l) are calculated relative 
to end I. Thus, the matrix equation (A-l) can be rewritten for = U. = 0, 

that is, 


'F7 > 


"AE/£ 

0 

O 

O 


'U7 ' 

Jf8 

y- 

0 

12El/£ 3 

O 

-6El/£ 2 

< 

U8 1 

|T10 


0 

0 

JG/£ 

O 


U10| 

T12> 


_ 0 

-6El/£ 2 

O 

4El/£ 


.U12> 
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and 


'FI' 


r F7 ' 


° 1 

If2 


F8 

> + < 

o 1 

T4 

T10 

1 

0 

^T6/ 


^T12> 


iF2 • t, 


U10 

V 



(a) Beam forces. 


T 10 



(b) Beam deflections. 

Figure A-l.- Beam deflections and beam forces. 
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APPENDIX B 


COMPUTER PROGRAM 

The mathematical model developed in the previous chapters was incorporated into 
a computer program so that the staging maneuver could be simulated as a function oi 
time. The program was written in FORTRAN V language and was run on the 
UNIVAC 1108 digital computer. The numerical integration routine used in the program 
was a variable step size Adams-Moulton method utilizing a fourth-order Runge-Kutta 
procedure for the required starting values. A description of the integration procedure 

is presented in reference B-l. 
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